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In this paper we present the main results of our investigation of the cqqq single- 
charm scalar tetraquarks and their SU(3)f representations: 165, 85, 6^ and 3^. We 
use the Fermi-Breit interaction Hamiltonian with SU(3) flavor symmetry breaking to 
determine the masses of the single-charm tetraquarks. We also discuss mass spectra 
obtained from meson and baryon mass fits. The mass spectra are very similar to those 
obtained with Glozman-Riska hyperfine interaction, and they indicate that some of the 
experimentally detected states may have tetraquark nature. 

^ I Key words: nonrelativistic quark model, hadron mass models and calculations, 

I light quarks, charmed quarks. 

«D . PACS: 12.39.Jh, 12.40.Yx, 14.65.Bt, 14.65.Dw 



1. INTRODUCTION 

The possible existence of four-quark states for light flavor dimensions, as well 
as some predictions for tetraquark spectroscopy, was first suggested by Jaffe H]. In 
Ref. f2| it is also provided a framework for a quark-model classification of the many 
two-quark-two-antiquark states. In Ref. |i3j| the energies of diquonia q^{(f') with or- 
bital angular momentum L = are calculated and compared to the treshold energies. 
The first observations of the scalar charmed mesons have been reported [?-'^. In 
Refs. |[9UTn some observed mesons (as for instance D+(2317)) are explained as a 
scalar cs systems. Many studies appeared in the past on D+(2317) ( |[T2] - [T4l and ref- 
erences therein). Possible interpretations include tetraquark states or molecular DK 
or DK* states |[T5l . The origin of the lightest scalar mesons, in particular, the light 
cr-meson is given in |[T6l in the framework of the instanton liquid model of the QCD 
vacuum. Terasaki and Hayashigaki II17II18II have investigated the decay rates of the 
members of the same multiplet in class of four-quark mesons. Also, in Refs. |[T9] - l2T]| 
it is shown that the existence of some exotic states with the four-quark structures 
might be expected. There are more results indicating that a diquak-antidiquark struc- 
ture is acceptable for some observed states Ii22ti26ll . In Refs. Il27ll28l a mixture of 
conventional quark-antiquark states and four-quark components is considered. Also, 
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there are predictions for some hidden charm states to be tetraquarks. For example, 
in 1291 the masses of the ground state heavy tetraquarks are calculated in the frame- 
work of the relativistic quark model, based on the quasipotential approach in quan- 
tum chromodynamics. These authors found that some exotic meson candidates can 
be tetraquark states with hidden charm. A concise overview of mesons with heavy 
quarks including charmed mesons and charmonium (or charmonium-like) states is 
given in 1301 . The mass spectrum of the scalar" hidden charm and bottom tetraquark 
states is studied in |[3TI . 

In this paper we perform a schematic study (two quark interaction) of the 
masses of the single-charm cqqq tetraquarks in the SU(3) flavor representations. We 
consider states with the spin-parity quantum numbers = 0"*" (J = total angular 
momentum, P = parity). Using the colored version of the Fermi-Breit (FB) hyper- 
fine interaction (HFI) Il3l [32] - l35l we investigate the possible four-quai"k structure of 
these mesons. According to the mass constraint, experimentally indicated states can 
contain only one c-quark. If each of them contains no heavy quarks, mass is too low, 
and if it contains one heavy quark (b or t) or two c-quarks mass will be too high. 
Only two possible solutions are: cqqq that we investigate and cqqq that is analyzed 
by Liu et al. |[23l . 

We showed |[36l that the constituent quark masses are very sensitive to the 
system in which they are contained, and their values differ less or more in differ- 
ent systems. That is why we analyzed tetraquark masses using independent meson 
and baryon fits for constituent quark masses. These fits satisfy Feynman-Hellmann 
theorem for FB interaction |[37l[38l . Constituent quark masses in tetraquarks are 
somewhat different for both cases: quark masses from meson fit and from baryon 

We choose relatively simple FB HFI because it nicely satisfies Roncaglia in- 
equalities for mass differences II37II38II . We deal with the open charm states. Since 
our intention is to investigate the possible tetraquark nature of 27 cqqq states, it is not 
necessary to develop an advanced relativistic quark model. That is why masses of 
the scalar charmed tetraquarks are discussed in the framework of the nonrelativistic 
quaik model, in which the mass of a hadron is considered to be the sum of the con- 
stituent quark masses and contributions of the FB HFI. As it will be shown later, our 
simple model is quite sufficient for investigating wave functions and masses of these 
states. Also, this model is used to get estimates of the theoretical meson and baryon 
masses. We prefer to deal with four quarks instead of diquarks because our system 
under consideration consists of one heavy and three light quarks. 

In this paper we present detailed calculation of masses using FB interaction, 
and it is the first time this interaction is applied to 27 states of cqqq tetraquarks in 
order to derive formulas for masses. We also compai^e these results with another 
phenomenological interaction: Glozman-Riska hyperfine interaction (GR HFI) BOl . 
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2. ANALYSIS AND METHOD 



We discussed single-charm tetraquarks in Ref. BTI . where the flavor wave 
functions and masses of cqqq tetraquarks ai^e calculated using GR HFI. Now we 
present detailed calculation of masses of cqqq tetraquai^ks using FB HFI. 

We analyze the tetraquark states with one charm quai^k which is singlet under 
the transformation of SU(3)f- There are four multiplets according to product: 3® 
3(8)3(8)1 = 155' + 85 + 6yi + 3a, i.e. two anti-triplets, one anti- 15-plet and one sextet. 
Young diagrams for these SU(3)f multiplets, as well as the weight diagrams, can be 
found in Figs. 1 - 5 of flU . In these tetraquark multiplets, all states have a chaim 
number equal to 1. Labels for all 27 states are the same as in Refs. 1141 [14211 . These 
labels are taken only by analogy with baryons, but of course they are not baryons. 
Many authors use the same labeling |[T9ll20l . The strong FB HFI Hamiltonian 131 
may be written in the following form: 

H^B = C^^{Xfxf), (1) 

i>j=i 1 

were ai are the Pauli spin matrices, Ap are the color Gell-Mann matrices and C is a 
constant proportional to strong hyperfine structure constant Oc. Hamiltonian ([T]) has 
explicit color and spin exchange dependence and implicit (by way of quark masses) 
flavor dependence. Its contribution to tetraquark masses is: 

m^^FB = W\{x\Hfb\x)W) , (2) 
where x denotes the spin wave function and u - flavor wave function. For total masses 
nil, we have: 

rrii, = mu,o + 'mu,FB, (3) 
where m^fi are masses without influence of FB HFI. 

Here we have to mention that the mixing of states is taken into account. In 
Table I of the paper BTI the four-quark content and quantum numbers of scalar- cqqq 
tetraquarks are given. One can see from that table which states mix due to the same 
quantum numbers. So we took into account the mixing when calculating masses of 
these tetiaquai^ks. 



3. RESULTS AND DISCUSSION 



In our model (total spin = 0), the corresponding symmetric xs and antisym- 
metric XA spin functions have the following forms: 

I X5) = ^ \nn + im + tut + tttt -2 tttt -2 tttt) , (4) 

IXA) = ^|ttti-ttti-tUt + itit>. (5) 
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The spin wave functions are symmetric (xs) or antisymmetric (xa) under in- 
terchange between the pair of quarks and between the pair of antiquarks (not for 
quaitantiquark interchange). Tetraquarks ai^e bosons, i.e. they have integer spin (for 
scalar tetraquarks it equals 0). For the open charm system, cqqq, the interaction be- 
tween the light quarks q and the c quark is suppressed in the heavy quark limit B3l . 
Thus, in the first approximation, three light quarks are decoupled from the heavy 
quai^k and there can be considered states compounded by u, d and s as triquarks or 
color nonsinglet baryons (in the bound state) B3l . Symmetry of their total wave 
functions is determined as in the case of fermions because there are three quarks in 
SU(3) group: qqq, where q = u, d, s. In the case of fermions the total wave function 
has to be antisymmetric, as well as the color state function in the case of all hadrons, 
and therefore the particles from multiplets 15s and 3s have the symmetric spin and 
flavor wave functions (xs^i^s)^ while the particles from multiplets 6^1 and 3a have 
the antisymmetric spin and flavor wave functions (xa, i^a)- 

The calculation of FB contribution to tetraquark masses will be described in 
more details. First we will explain how to calculate the ffiUj products of Pauli spin 
matrices for each pair in the scalar system of four quarks. These products are the 
expected values of spin matrix elements. We use the spin operator eigenvalues for 
triplet and singlet states: (-1-1) and (-3), respectively. We also use the following values 
of the matrix elements (Jia2: 

(ni<?i^2iti) = -i = ati'?i'?2iit) 

(ttl^l^2|tt) = l = (Ul^l<?2|U> (6) 

ati^i<?2iti) = 2 = (ni<5^i'5^2iit). 

In addition to that, for total spin S and its projection mg, for symmetric state it 
holds: 

\S = l,ms = l) = \n) 
|5 = l,m, = 0) = ^|n + it) (7) 
\S = l,ms = -1) = III) , 

from where we get + ||) = = l,ms = 0). 

For antisymmetric state it holds: 

|5 = 0,m, = 0) = -^|n-it>, (8) 
and combining expressions ^ and ©, we get: HJ,) = ^ {\S = 1) + \ S = 0)) and 
\lt) = ^i\S = l)-\S = 0)). 

Then we calculate symmetric {xs\^i^j\xs) and antisymmetric matrix ele- 
ments {xa\ ^i^j \xa) for the following pairs: qiq2, qm, qm, 92^3, ^2^4, Mi, and 
apply them to relations for spin wave functions dH) and For example, for qiq2 
pair, we use only spins of the first and second particle, therefore the other two spins 
e.g. of the third and forth particle are used only to determine which addends are 
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non-zero. In case of symmetric spin wave function it leads from eq. Q to eq. (ITOl ): 

{Xs\^1^2\Xs) = 

= ^ ((tll2 + Ilt2) t3i4 + (tli2 + ilt2) i3t4 + 

-2tlt2i3i4-2ili2t3t4|x (9) 
X {B1B2) I (tli2 + ilt2) + (tli2 + ilt2) i3t4 + 

-2 tlt2i3i4 -2 Ili2t3t4) 



(X5|'7l(T2|XS> 

(10) 



4^ ( ( (ti + it) + (ti + it) -2 tt -2 ii I X 



X {0102) I (ti + it) + (ti + it) - 2 tt -2 ii)) . 

Nevertheless, we see that addends labeled with I and II in eq. (ITOl ) cannot be com- 
bined into non-zero matrix elements because the spins of the third and fourth particles 
are t| in case I, and \;\ in case II. So we derive: 

{Xs\dlS2\Xs) = 

^ {y/2{S = l,m, = 0| (<Tia2) \S = = 0) \/2+ 
+V2{S = l,m, = 0| (ai<T2) |5 = l,m, = 0) ^/2+ (11) 
+ (-2) {S = l,m, = 1| (ai^a) |5 = l,m, = 1) (-2) + 
+ (-2) {S = l,m, = -1| (ai<72) \S = l,m, = -1) (-2)) , 
from where we obtain this result: {xs\ <?i'?2 \xs) = ]^(2-l + 2- l + 4- l + 4-l) = l. 
In this way we get the results for all quai^k pairs qiqj: 
{xs\ 5\a2 \xs) = 1 = {xs\Szd^ \xs) 

kXs\ o'l<?3 \XS) = -2 = (X5| '?2'?4 \XS) |-j2) 
{xa\ S\a2 \XA) = -3 = (xaI <?3'?4 Ix^) 
(X^l '?i'?3 \XA) = = (xaI 520^ \XA) ■ 

Table 1. 

The products (Gell-Mann matrices for color SU(3)c) and the products BiBj (Pauli spin 
matrices) for symmetric and antisymmetric multiplets. 



multiplet 










QiQj 






15s and 85 


A1A2 = A3A4 = 


8 
3 


A1A3 — 


A1A4 = 


A2A3 = 


A2A4 = 


4 
3 




(71IT2 = <?3'?4 = 


: 1 


aia-s = 


(71(74 = 


^ cr2(73 = 


O20-4 = 


-2 


6a and 3^ 


A1A2 = A3A4 = 


8 
3 


AiA3 = 


A1A4 = 


A2A3 = 


A2A4 = 


4 
3 




(T1CT2 = (?3i?4 = 


-3 


(71(73 = 


= (71(74 = 


= 0-20-3 = 


= 0204 = 


:0 



In Table [T] we give the products of Gell-Mann matrices for color SU(3)c and 
the products of Pauli spin matrices. When we put values for ffiffj and A/Aj from 
Table [TJinto equation we get the following FB HFI contributions for symmetric 
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Table 2. 

Masses of scalar cqqq tetraquarks distributed in SU(3)f multiplets, with mixing between states with 
the same quantum numbers. m^^Q are tetraquark masses without influence of FB HFI and m^^ pB are 
FB HFI contributions to tetraquark masses. 

multiplet cqqq m^fi (m„ = m^) m,,,FB = "^d) 



3 \ msTric rriums muiric mi. 



Ha IrUu+ms + rUc |C ( + ; 



1 



3 V rnsTUc 

A 3mu + rric |C f ^ H ^ 

S 2m„ + ms + mc |C(— 4;- + 



* c 3 V rriumc ~ mums m.sm.c mf, 



155-35 2m„ + m, + m,;3m, + me |C + ; |C + 



mf. m„mc / ' 3 \ m^mc m- 



64 So 2m„ + m, + m^. 8C ( — ^ 1 ^ — | 

"A * use ymums mumc J 



Q Irriu + nis + rric 8C I + 



1 I 1 



3a Ds 2m„ + m, + mc 8C — — + 



mf. msmc 



1 I 1 



mums mumc 



Ga-3a D 3m„ + m^; m„ + 2m, + me ( -1- + -i- ; ^ + ^ 



mums msmc I ' \mf' mum, 



and antisymmetric multiplets: 

8^/ ,1 1 , V 

ruy^FB^s = -^C{vs\ \ Ws) + 

6 mim2 m^m^ 

+lc{i^s\^— + ^— + ^— + ^—Ws), (13) 
6 niim^ mim^ 11121713 111,21714 

for symmetric multiplets 15^ and 35, and 

mu,FB,A = 8C{l^A\ \ Wa), (14) 

11111112 in^m^ 

for antisymmetric multiplets 6a and 3a- The flavor wave functions vs and ua for the 
scalar cqqq tetraquarks are given in Table II of Ref. ATI . The masses of tetraquarks 
predicted from our model are given in Table |2] The mixing between states with the 
same quantum numbers is included. There is mixing between states from symmetric 
multiplets 15s and 3s, and also between antisymmetric multiplets 6a and 3a, while 
symmetric and antisymmetric multiplets do not mix with each other. We first show 
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our predictions for spectra when FB HFI is not included (the third column of the 
table) and then we show the FB HFI influence (the forth column). 

Here we study mass spectra of single-charm tetraquarks using FB HFI in sche- 
matic approximation (two-paiticle interaction). The masses of constituent quarks 
(= nid), rus and rric and the constant are calculated from fitting of the equations 
for meson and for baryon masses, with FB interaction included, to the experimental 
meson and baiyon masses ll36l . As one can see from Table 2 in Ref. [44.1 or Table 2 in 
Ref. P31 and from references therein, our predictions for constituent quark masses 
are similar to masses obtained using different phenomenological models. We use 
masses for these mesons: light pseudoscalar mesons vr, K, light vector mesons p, 
K* , w, 93, charmed mesons D, D* and strange charmed mesons Ds, D^. We did not 
calculate 7] and rj' contribution because of their mixing and because they cannot be 
described within such a model. Also, we use masses for these baryons: light baryon 
octet N, S, E, A, light baryon decuplet A, S*, S*, $7 and heavy baryons Sc, Ac, S*, 
f^c- For each set of equations, the minimized values for masses are calculated by 
formula (14) given in Ref. ||4T1 . The corresponding experimental masses are taken 
from the "Particle Data Group" site: http://pdg.lbl.gov Il46l . 

From the fit of all meson masses (when mesons with two c-quarks are ex- 
cluded) we obtained the following values |[36l (Table III): ruu = rrid = 314.75 MeV, 
ms = 466.80 MeV, = 1627.31 MeV and the constant C"" = 1.5546 x lO^MeV^. 
These values for quark masses here we use for calculating the tetraquai^k masses. 
For the constant C we use value C^^^"^"" which is different from (7""=*°". We fit- 
ted (7*'^*''" in that way to obtain the mass of the lowest state from 3a equal as the 
D+(2317) meson: C*^*™ = -5.80 x 10*^MeV^. From the fit of all baryon masses 
(when barions with two c-quarks are excluded), in |[36i (Table III) we obtained: m„ 
= md = 365.69 MeV, = 530.08 MeV, rric = 1700.17 MeV and the constant = 
1.2513 X lO'^MeV^. Here, we calculated C*^*™ = -11.90 x lO^MeV^ (in that way 
to obtain the mass of the lowest state from 3a equal as the D+(2317) meson). 

The values of theoretical masses of some mesons and baryons with FB HFI 
included, using constituent quark masses and the constants C™, given in Table 
III of Ref. |[36ll (the upper rows which correspond to FB HFI) we present in Tables [3] 
andlH If we compare Tables [3] and |4] we can notice that theoretically obtained masses 
of baryons using FB HFI are in better agreement with the corresponding experimental 
masses than masses of mesons. 

The values of masses (in MeV) of scalar cqqq tetraquarks, obtained from the 
meson fit, are given in Table|5l and their masses obtained from the baryon fit are given 
in Table (6] These tables with tetraquark masses show us where the 27 tetraquark 
states are expected. There are uncertainties in calculating masses because we use the 
model with schematic interaction. 

The tetraquark mass spectra are given in Figs. [T] and [2] (labels are the same as 
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Table 3. 

The values of theoretical masses mm (in MeV) of some mesons with FB HFI included, when we use 
constituent quark masses and the constant C™ obtained from the meson fit (Table III of Ref. 1361 ). 
rUm^ are experimental masses 1461 and Amm is the absolute difference between these two values. 



meson 


mm (MeV) 


m^'P (MeV) 


Arrim (MeV) 


vr 


159 


140 


19 


K 


464 


494 


30 


P 


786 


776 


10 


K* 


887 


892 


5 




786 


783 


3 




1005 


1020 


15 


D 


1851 


1869 


18 


D* 


1972 


2010 


38 


Ds 


2033 


1968 


65 




2015 


2012 


3 



Table 4. 

The values of theoretical masses mt, (in MeV) of some baryons with FB HFI included, when we use 
constituent quark masses and the constant obtained from the baryon fit (Table III of Ref. 1361 ). 
mf^^ are experimental masses 1461 and Am^ is the absolute difference between these two values. 



baryon 


rub (MeV) 


ml''^ (MeV) 


Amb (MeV) 


N 


957 


940 


17 


S 


1179 


1190 


11 




1319 


1315 


4 


A 


1121 


1116 


5 


A 


1237 


1232 


5 


S* 


1373 


1385 


12 


"* 


1513 


1530 


17 


n 


1657 


1672 


15 


So 


2438 


2455 


17 


Ac 


2291 


2285 


6 



in Ref. HTTl ). It can be noticed that HFI determines mass splitting in the spectrum, 
i.e. its fine structure. And for mixing it can be said that it separates two states. When 
comparing Fig. [T] (or Fig. ^ from this paper with Fig. 6 from BTI and Fig. 1 
from 1421 . it can be noticed that GR HFI reduces masses (except for IS^ - 85 mixed 
states) more than FB HFI, but that difference is not so significant. One can see that 
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Table 5. 

The values of masses (in MeV) of scalar cqqq tetraquarks distributed in SU(3)f multiplets, with 
mixing between states with the same quantum numbers, obtained from the meson fit. q (MeV) are 
tetraquark masses without influence of FB HFI, pB (MeV) are FB HFI contributions to tetraquark 
masses and mj/ (MeV) are the total tetraquark masses. 



multiplet 


tetraquark 


m^,o (MeV) 


mi.,FB (MeV) 


m,y (MeV) 


15s 




2876 


-150 


2726 






2724 


-176 


2547 




A 


2572 


-186 


2385 




S 


2724 


-94 


2629 


T5s-3s 


D,(T5s - 3s) 


2724; 3028 


-176; -91 


2547; 2936 




D(15s-3s) 


2572; 2876 


-186; -59 


2385; 2817 


6^ 




2724 


-406 


2317 







2724 


-529 


2194 


3a 




2724 


-406 


2317 


Qa-3a 


D(6a - 3a) 


2572; 2876 


-377; -559 


2195; 2317 



3.4-r 

3.2- 




2.0-1 1 1 1 1 1 

15s -^S -^A 

Fig. 1 - Tetraquark mass spectrum from the meson fit without (left column) and with (right column) 
FB HFI, both with SU(3)f symmetry breaking. 

FB HFI reduces the obtained masses and causes splitting between S and T^s (Ds) 
in 15s and between E^j and in 6a, and GR HFI only causes splitting between 
and Q. Also, there is difference for 15s — 3s mixed states when comparing with GR 
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Table 6. 

The values of masses (in MeV) of scalar cqqq tetraquarks distributed in SU(3)f multiplets, with 
mixing between states with the same quantum numbers, obtained from the baryon fit. q (MeV) are 
tetraquark masses without influence of FB HFI, m,y pB (MeV) are FB HFI contributions to tetraquark 





masses and tcvu 


(MeV) are the total tetraquark masses. 




multiplet 


tetraquark 


my^Q (MeV) 


m^_FB (MeV) 


rriy (MeV) 






3126 


-234 


2892 






2962 


-273 


2689 




A 


2797 


-288 


2509 




S 


2962 


-157 


2805 


T5s-35 


D,(T5s - 3s) 


2962; 3290 


-273; -148 


2689; 3142 




D(15s - 3s) 


2797; 3126 


-288; -102 


2509; 3024 


6a 




2962 


-644 


2317 






2962 


-818 


2144 


3a 




2962 


-644 


2317 


6a -3a 


D(6a - 3a) 


2797; 3126 


-597; -865 


2200; 2261 



3.4 




2.0 J + f \ + — 

15s -^S -^A 

Fig. 2 - the same as Fig. 1, but for baryon fit. 



HFI, but the forms of tetraquark spectra with FB and GR interactions are similar. 
It is interesting to note that the spectra are similar- although the one interaction is 
color-spin, and the other one is flavor-spin. 

From spectrum it is possible to identify D+(2317) as the lowest state in multi- 
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plet 3 A- It agrees with the identification of D^(2317) with tetraquark state which was 
considered in Ref. |[47| . In Ref. BTl . for D+(2632) it was claimed to be a candidate 
for a tetraquai^k D^, but they proposed different total angular momentum (J^ = 2~^). 
In the present paper, the authors consider scalar tetraquark only. Also, from spectrum 
we can see that FB HFI reduces the obtained masses for all states. 

States D+(2632) H and D0(2308) ^ are mixed states: D+(2632) is from 
lbs - 3s and D°(2308) from 6a - 3a mixing. Because of that mixing, their flavor 
wave functions ai^e given only in a first approximation [H and therefore calculations 
are not sufficiently precise. Therefore these two states have theoretical predictions 
which are not the same as the experimental ones. But nevertheless, for meson fit, 
their experimental masses are still between the two values we obtained in Table [5] 
D+(2632) has the mass between 2547 and 2936 MeV, while D°(2308) has the mass 
between 2195 and 2317 MeV. Obtained theoretical masses for baryon fit in case 
of D°(2308) is lower than it it expected and for D^(2632) is approximately in the 
expected range. If we compare Tables |5] and [6] and Figures [T] and |2] we can conclude 
that FB interaction gives lower total masses. FB interaction that we applied needs 
some further improvement. Probably the biggest discrepancy between theoretically 
obtained and experimentally detected state is because of sensitivity of constituent 
quark masses to the system in which they are contained (meson, baryon, tetraquark). 

4. CONCLUSIONS 

We calculated the mass spectra of two-quark-two-antiquark system. We used 
the observed meson masses, taken from "Particle Data Group" B6l . to obtain con- 
stituent quark masses by way of mass fit. Applying the Hamiltonian ([T|l to the 
constituent quarks, we obtained the theoretical meson, baryon and single-charm 
tetraquark masses with FB contribution included. For the first time the FB HFI, 
which is color-spin interaction, is applied to 27 states of cqqq tetraquarks in order to 
derive formulas for masses and these are new results. Also, when comparing with 
GR HFI, which is flavor-spin interaction, it is interesting that we also obtained simi- 
lar results for tetraquark masses, like in the case of meson and baryon masses. In our 
constituent quark model, FB HFI is calibrated in mesons and baryons, and then used 
in tetraquarks. 

If the tetraquark states we have studied contained no heavy quarks, their mass 
would be too low; if they contained one heavy quai^k (6 or t) or two c-quarks, mass 
would be too high. According to that and to the mass constraint of the experimen- 
tally indicated states, these states can contain no more than one c-quark. We also 
showed |[36l that the constituent quark masses are very sensitive to the system in 
which they are contained, and their values differ less or more in different systems. 
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That is why we analyzed tetraquark masses using meson and baryon fits for con- 
stituent quark masses, independently. Constituent quark masses in tetraquarks are 
somewhat different for both fits. We choose relatively simple FB HFI because it 
nicely satisfy Roncaglia inequalities for mass differences II371I38II . 

Symmetric 15s-plet mixes with the 3s-plet ideally (Dg and D states). This 
mixing splits the two states into a heavy (hidden strangeness) and a light one. Also, 
antisymmetric 6y!i-plet mixes with 3yi-plet (lowest mass D is the ideal mixture of 
these antisymmetric multiplets, while the lowest is pure 3yi-plet). From the spec- 
tra of single-charm tetraquark masses with Fermi-Breit HFI, it can be noticed that 
this interaction implies no flavor dependent splitting among multiplets. 

Flavor wave functions of the mixed states are given only in a first approxima- 
tion (see Ref. H]). The mixing of the states also changes the properties and shifts 
masses from the theoretical predictions. For instance, possibly tetraquark states 
D'^(2308) and D^(2632) in our case are mixed states and the calculations of their 
masses are not sufficiently precise. D+(2632) appears as a mixed state from mixing 
of multiplets I5s and 3s, and 0*^(2308) would be from mixing 6a and 3^- Accord- 
ing to our results, all three states D+(2632), D°(2308), D+(2317) might have the 
tetraquark nature. We gave the contribution of FB HFI to cqqq tetraquark masses and 
also we calculated tetraquark masses and compared them to GR HFI. 

As it can be seen from Tables [5] and [6] and Figs. \T\ and |2l FB HFI reduces the 
obtained masses and causes splitting between S and in 15s and between and 
17 in 6^. Besides, the spectra obtained from different fits have a similar aiTangement 
of particles. Probably the biggest difference between theoretical and experimental 
states is due to: 

(i) sensitivity of constituent quark masses on systems in which they are contained, 

(ii) wave functions for the two detected experimental states 0^^(2308) and D+(2632) 
are calculated only in a first approximation, and obtained mases are not precise, 

(iii) FB HFI is not the completed HFI. 

More experimental searches for detection of other cqqq members are needed in 
the future. 
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